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Spin relaxation in the impurity band of a 2D semiconductor with spin-spUt spectrum in the ex- 
ternal magnetic field is considered. Several mechanisms of spin relaxation are shown to be relevant. 
The first one is attributed to phonon-assisted transitions between Zeeman sublevels of the ground 
state of an isolated impurity, while other mechanisms can be described in terms of spin precession in 
a random magnetic field during the electron motion over the impurity band. In the later case there 
are two contributions to the spin relaxation: the one given by optimal impurity configurations with 
the hop-waiting time inversely proportional to the external magnetic field and another one related 
to the electron motion on a large scale. The average spin relaxation rate is calculated. 

PACS numbers: 71.55.Jv, 71.70.Ej, 72.25.Rb, 85.75.-d 



Spin dynamics in semiconductors has attracted much 
attention in the last decades 1, In particular, a 
number of experimental 0, 0, 5, [a, [3| and theoretical 
@, i, M, [H works are devoted to the investigation of 
spin relaxation in the impurity band of a semiconductor. 
An increasing interest to this problem is motivated by 
experimental observation of up-to-300ns spin lifetimes in 
n-doped bulk GaAs and GaAs/AlGaAs heterostructures 
[1) |j) S @1 J which makes them good candidates for the 
use in possible spintronics applications. Yet, a consis- 
tent theory of spin relaxation in the impurity band is 
still to be developed. Depending on the donor concen- 
tration, spin relaxation in the impurity band might be 
driven either by hyperfine interaction or spin-orbit cou- 
pling. Since the nuclear spin relaxation time is typically 
very long, hyperfine interaction can be treated as a ran- 
dom-in-space static magnetic field with the associated 
spin precession frequency ujn ^ A/^/N, where A is the 
hyperfine coupling constant and N is the number of nu- 
clei within the volume occupied by the wave function [l^ 
(the directions of the random magnetic field for electrons 
located on different impurities are not correlated) . In the 
case of spin-orbit coupling, the associated spin precession 
frequency is a power function of the electron momen- 



tum p 
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(in the 2D case, Wp is linear in p). 



As a result, spin-orbit coupling leads to spin rotation in 
the process of phonon-assisted hops from one impurity 
to another by the angle </) « u)pgAr/vQ, where Ar is the 
distance between impurities and po = mvg is the under- 
the-barrier momentum. There are several mechanisms 
of spin relaxation in the impurity band. As in quantum 
dots, spin relaxation might be driven by phonon-assisted 
transitions between Zeeman sublevels of the ground state 
of separate impurities. Other mechanisms, that involve 
electron hops from one donor to another, are specific for 
the impurity band. In Ref. the spin relaxation rate 
was estimated as: 

1/ts ^ UjjjThc, 1/ts ^ Cj)^ /The (1) 

for the case of hyperfine interaction and spin-orbit cou- 



pling respectively and the characteristic hop waiting time 
The was assumed to depend only on the average distance 
between impurities. These equations are based on the 
classical picture of the angular spin diffusion in a random 
magnetic field (in the case of hyperfine interaction, the 
direction of spin precession changes randomly after each 
hop; in the case of spin-orbit coupling, the spin rotates 
by a certain angle in a random direction in the process 
of a hop). However, this approach does not account for 
the exponential variation of the hop waiting times: 



Thi = Toexp (2Ar/a) , 

Th2 = tq exp (2Ar/a + A£/T) 



(2) 
(3) 



for phonon emission and absorbtion respectively (here 
Ar is the distance between impurities, AS is the dis- 
tance between the energy levels, a — eh?/2me^ is the 
Bohr radius, and T is the temperature). The main con- 
sequence of such inhomogeneity is that it is impossible 
to introduce an universal time scale for the system under 
consideration. This fact is confirmed by about ten-fold 
decrease of the experimentally measured spin correlation 
time in the bulk GaAs at the crossover from hyperfine- 
interaction-induced to spin-orbit-induced spin relaxation 
(see Fig. 3 in Ref. 0). The effects of the inhomogene- 
ity on the spin dynamics in the absence of the external 
magnetic field were considered in Refs. [I,[l3| for the sys- 
tem with spin-split spectrum. In particular, it was found 
that there are two essentially different contributions to 
the spin relaxation: the one related to electron hops over 
the pairs of impurities with the size of the order of the 
Bohr radius and another one related to the motion over 
a long distance. 

In this letter, we calculate the average spin relaxation 
rate for the mechanisms discussed above in the presence 
of the external magnetic field. The use of the averaged 
relaxation rate is justified if the relaxation is slow enough 
so that an electron can walk over a large distance dur- 
ing the spin relaxation time (in the opposite case the 
spin relaxation is governed by escape from the regions 
with slow relaxation to the regions with fast relaxation 
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[lot). The corresponding condition is ts ^ re (here 
Tc = Toexp(C^o) is the hop waiting time for so-called 
critical bond, = \/^L'^W/a'^T, C is the dimensionless 
coefficient [131, ^ — e^/e-^d is the characteristic width 
of the impurity band, and = is the average 

distance between impurities). We assume that spin pre- 
cession in the external magnetic fields is sufficiently fast 
^o'l'S ^ 1 (here fio is the spin precession frequency in 
the external magnetic field B). In this case, the compo- 
nents of the spin perpendicular to the magnetic field are 
suppressed due to fast precession, and hereafter they will 
be neglected. We assume that the temperature is suffi- 
ciently small T <C W, so that we can neglect activation 
to the conduction band. We also assume that hflo <C T, 
neglect electron-electron interaction, and treat donors as 
2D coulomb centers. 

Our point is that over a wide range of magnetic fields 
the relevant time scale The for the problem under consid- 
eration is given by: 



(4) 



Indeed, a common feature of the relaxation mechanisms 
based on the angular spin diffusion in a random magnetic 
field is that they are suppressed by applying a longitu- 
dinal magnetic field with the associated spin precession 
frequency larger than the inverse correlation time of the 
random magnetic field. In the simplest case of a pair of 
impurities with the hop waiting times r/ii = t;j2 = 
(AS <^ T), the spin relaxation rate is proportional to 
AfpTfi/ (l + ripT^), where Ail is the spin precession fre- 
quency in the random magnetic field (in the case of hy- 
perfine interaction AO « lon] in the case of spin-orbit 
coupling AQ w i^o^j as shown below). The contribution 
of the pairs to the spin relaxation increases exponentially 
with Ar for < I/Hq and decrease for Th > 1/rio- Tak- 
ing into account Eqs. H]) and we can estimate the 
spin relaxation rate on the pairs of impurities as: 



(5) 



for the case of hyperfine interaction and spin-orbit cou- 
pling respectively (here z/ ^ {a/ Ldf T/W is the share of 
the optimal pairs) . We also consider other relevant mech- 
anisms of spin relaxation that are related to the electron 
motion over a long distance and spin-flip processes. 

Let us proceed to the rigorous formulation of the prob- 
lem. First, we consider the system with spin-split spec- 
trum. The Hamiltonian of the system is 



(6) 



where Hq and ifph are the Hamiltonians of an electron 
and phonons respectively, and the last term on the right- 
hand side describes electron-phonon interaction. The Ha- 
miltonian of an electron is given by: 



2m 



U{r) + h<Tno/2 + h(Tap/2mLs, (7) 



where U (r) is the impurity potential, Ls is the length 
characterizing the strength of the spin-orbit coupling, a 
is the dimensionless tensor with the components of the 
order of unity, and cr is the vector of Pauli matrices. The 
last term on the right-hand side is a combination of the 
Bychkov-Rashba spin-orbit coupling [l^ and Dresselhaus 
spin-orbit coupling averaged over the electron motion in 
the direction perpendicular to the quantum well 3, l3] ■ 
The Hamiltonians of phonons and electron-phonon inter- 
action are given by: 



q 



(8) 
(9) 



where 6q and bq arc phonon creation and annihilation 
operators, q is the phonon wave vector, s is the sound ve- 
locity, Cn is the coefficient characterizing the strength of 
the electron-phonon interaction, V is the system volume, 
and n = ±1 for deformation and piezoelectric phonons 
respectively. For the following consideration it is conve- 
nient to make a transformation, which cancels spin-orbit 
coupling to the first order in parameters l/Ls and Hq: 



fjl _ giCTQr/2Ls^g-iCTar/2Ls 



(10) 



As a result. 



HL^ ^ + U{y) + nno(T/2 + n [Oo X ar/Ls] (t/2, (11) 
2m 

while i/ph and i/c-ph are not modified. 

First, let us consider the spin relaxation caused by pho- 
non-assisted transitions between Zeeman sublevels of the 
ground state of an isolated impurity. In this case the 
spin relaxation rate coincides with the transition rate ob- 
tained by the Fermi golden rule: 

— = ^ / iVqWq S mo - hsq) Vdq/{2Tr)\ (12) 
TS n J 

where Wq — 1(5'+ |exp (iqr)| \t'_)|^, A^q is the phonon 
occupation number, 'Si = \I>o± + 55'±, 



(13) 



are the eigenfunctions of an electron at a Coulomb center. 



(5*± =(£o-Hc) [ar/Ls x nno/2] ct^o± (14) 



are the corrections due to the last term on the right-hand 
side in Eq. (fTTj) . |±) are the spinors, He = p^/2to— e^/er 
is the Hamiltonian of an electron at a Coulomb center, 
and £o — h^/2ma^ is the binding energy. In the case of 
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small magnetic fields ftoa <C s and fi-flo ^ T, we get: 
a aC, 



1 

TS 



T 

h \Ls So 



2 / X 5+71 
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(eo) = ^ |eo X ae^ 



(15) 
(16) 



where /q = (^'o| r (^£0 + ^c) r l^'o) ^o/a^ ~ 1 and eo 
is the unit vector in the direction of the external magnetic 
field. This dependence of the spin relaxation rate on the 
external magnetic field is in accordance with the results 
obtained in Ref. [l^ for the circular quantum dot with 
piezoelectric phonons {n — 1). 

Next, let us consider the spin relaxation on a pair of 
impurities caused by spin precession in the random mag- 
netic field. From the Hamiltonian (fTTjl one can derive an 
equation, describing spin dynamics: 



dS/dt=[{nQ + An{t)) X S]. 



(17) 



where AJl (t) = [fJo x ctr (t) /Ls] and the position of the 
electron r (t) takes two value: ri or r2 (here ri^2 are the 
positions of the impurities). To find the random mag- 
netic field correlator k {t) = {Aft (t) Aft (0)), we use the 
kinetic equation for an electron on a pair of impurities: 



drii/dt = -dn2/dt = n2/Th2 - ni/rhi 



(18) 



where ni^2 are the probabilities to find an electron at im- 
purity 1 and 2 respectively. Using the Green function of 
the kinetic equation, we get: 



n{t) 



Aft^ 



4cosr {AE/2T) 



exp {-t/Th) 



(19) 



where Aft = [TIq x aAr/Lg], Ar — Vi — r2, and l/r^i = 
^/Tfii + l/T'/i2- Treating the term proportional to Aft (t) 
in Eq. ([T7]) as a perturbation and using Eq. (fTOl) . wc get 
the following evolution equation for the component of the 
spin parallel to the external magnetic field: 



dSJdt- 



K{t')cos{not')Sii{t-t')dt', (20) 



The spin relaxation rate on a pair of impurities is 



l/rs(Ar,A£) 



An^ 



4cosr {AE/2T) 1 



(21) 



Depending on the strength of the external magnetic 
field, several regimes of spin relaxation can be realized. 
In the case fio < I/tq, the main contribution to the spin 
relaxation rate comes from the pairs of impurities with 
\A£\ < T and Tft sa l/VL^ (as follows from Eq. dUl), this 
contribution is proportional to the first power of the ex- 
ternal magnetic field). The average spin relaxation rate 
is 



l/rs = j l/Ts{AY,AE)dArdAElWLl. 



(22) 



Substituting Eq. into Eq. ([221), we get: 

2 



1 



64 



a a 



g(eo). (23) 



In the case ^Iq > I/tq, we can neglect the unity in the 
denominator of Eq. (I21|) . As a result. 



l/rs(Ar,Af) 



1 [eoxaAr/Ls] exp(— 2Ar/a) 



To 



l + exp{A£/T) 



(24) 

In this case the main contribution to the spin relaxation 
rate comes from the pairs of impurities with Ar < a and 
\A£\ < T. Substituting Eq. dHD into Eq. we get: 



1 __ 37rln2 I T / a a 
TS 4 To W lis Ld 



9 (eo) 



(25) 



Thus, the spin relaxation rate is saturated at large Qq- 

In deriving Eq. we assumed that optimal pairs 

are separated from the rest of the system, i.e. that an 
electron makes many hops over the pair before it leaves 
it. For this assumption to be valid, it is required that 
^qTc ^ 1. In the opposite case, an electron motion on 
a large scale gives the leading contribution to the spin 
relaxation rate. This contribution is proportional to the 
electron diffusion coefficient [l3| : 



1/ts ^D^ 1/tc. 



(26) 



The influence of the external magnetic field on the diffu- 
sion coefficient is well known ilT] . At low fields Rc ^ aS^o 
(here Rc is the cyclotron radius and a^o is the optimal 
hopping length [l7|), it can be described in terms of bend- 
ing of the tunnelling electron trajectory by the external 
magnetic field perpendicular to the quantum well, which 
effectively increases the distance between impurities: 



a^o {H) = a£,o 

As a result, 

1/ts {H) = (I/T5) exp 



1 + lf^^ 
60 \Rc 



'60 (i?c) ^° 



(27) 



(28) 



Next, let us briefly consider the spin relaxation caused 
by hyperfine interaction. In this case, Eq. pT|) still can be 
used with the replacement AD, — > Afl' = 2^^^ll}n. Using 
Eq. (|2ip and following the same procedure as before, we 
get: 



1 

TS 
1 
TS 



1 



= 7rln2 



T 

NVlIto W \L, 



LdJ \DoTo 



a 



(29) 
(30) 
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for the case fio < I/t'o and f^o > l/^o respectively. 

In the case of small magnetic fields ^IqTc <C 1 the spin 
relaxation rate can be estimated as 



l/rs = / Ah/NdP{t) 



(31) 



where P (t) is the probability for an electron to spend 
time t at an impurity. This probability can be replaced 
with the probability that an impurity is separated from 
the rest of the system by the length aln(t/ro) in the 
coordinate space and Tin (t/ro) in the energy space. As- 
suming that the form of the surrounding empty area is 
given by \Ar/a + A£/T\ < In [t/ro], we get: 



P(i) =cxp(-ln3 [i/ro] 87r/3^3) 



(32) 



Substituting Eq. ((32)) into Eq. (|3T|) and using saddle- 
point approximation, we get: 



1/rs « {7r/32f' (A^N) Co'/^exp (|v^) " (^3) 



The influence of the external magnetic field on the orbital 
motion can be accounted for by the replacement: 



(34) 



Here a^y£^/8^T is the size in space of the empty area at 
the saddle point. As a result, 



1/ts{H) = (1/Ts)exp 



60 (i?^) ^^"^ 



,3/2 



(35) 



To conclude, the theory of spin relaxation in the impu- 
rity band of a 2D semiconductor in the external magnetic 
field is presented. It is shown that spin precession in the 
external magnetic field enhances spin-orbit-induced and 
suppresses hyperfine-interaction-induccd spin relaxation. 
For spin orbit coupling, the relaxation rate is linear in B 
over a wide range of parameters, while the dependence 
on the direction of the external magnetic field is the same 
for all spin-orbit-induced relaxation mechanisms. For hy- 
perfinc interaction, the spin relaxation rate is inversely 
proportional to the external magnetic field. 
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